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1. Introduction 

The Pythagorean triples of integers satisfying x 2 + y 2 = z 2 have been 
studied and enumerated since Babylonian times. Since Diophantus, it has 
been known that this set of triples is related to the standard rational param- 
eterization of the unit circle, (fsrT) £rrj)- The Pythagorean triple solutions, 
which are relatively prime, have the elementary and beautiful characteriza- 
tion as integers x = m 2 — n 2 ,y = 2mn, z = m 2 + n 2 for relatively prime in- 
tegers to, n. One can also realize that the Pythagorean triples are related to 
the Gaussian integers, Z[i], the lattice in the complex numbers with integer 
coordinates, u = x + iy, where x, y, are integers. The Pythagorean equation 
N(u) = uu = (x + iy)(x — iy) = z 2 is an equation among Gaussians. It is 
perhaps, not surprising, that the Pythagorean triples are just the squares of 
the set of Gaussian integers, that is v = (to + ni) 2 = (to 2 — n 2 ) + (2mn)i 
gives the endpoint of the hypotenuse of a right triangle with integer length 
sides. This Gaussian square v obviously has a norm, N(v), which is a square. 
Both of these expressions secretly involve the double angle formulae for sine 
and cosine since the stereographic projection formula uses the central and 
chordal angles. 

By considering the action of the modular group, T = PSL^Z), by con- 
jugation on the set of all integer matrices M2(Z), we can blend together 
these two perspectives. We shall show (Corollary |3.2|) that the Pythagorean 
triples can be identified with an orbit of a subgroup (of index 6), T(2), 

which is generated by the images of U 2 = (I fj and L 2 = (l °X More- 

over T(2) is the free product of the subgroups generated by U 2 and L 2 



(Proposition |3.1|) . The consequence of the free group structure for T(2) is 
that the Pythagorean triples can be enumerated by words, for each k > 0, 
U ±2n kL ±2m kU ±2n k ^ 1 ... L ±2m lU ±2n lL 2m 0j for inte gers TO; / 0, n< + 0, 

< i < k. 

Using a breadth-first enumeration of these group elements we obtain the 
following curious result (see Figure HJ) . 

Theorem 1.1. The set of primitive Pythagorean triples has the structure 
of a complete infinite rooted ternary tree. 
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Figure 1 . Tree of Pythagorus 



2. Conjugation Action and Similarity 

Proposition 2.1. Any integer matrix X , satisfying X 2 = is X = ^ ^ 

for integers x, y, z satisfying yz + x 2 = and conversely. 

Proof. Suppose the integer matrix X, satisfies X 2 = 0. Since the matrix X 
also satisfies its characteristic polynomial, it follows that tr{X)X = det{X)I. 
Since it is impossible for a non-zero X to be a non-zero multiple of the 
identity matrix, it follows that tr{X) = det{X) = 0, and the conclusion of 
the lemma is immediate. ■ 

The group GL2(Z) of invertible linear transformations of Z 2 acts via con- 
jugation on all two by two matrices, M2(Z), by X — » TXT~ l . This action 
preserves the trace and determinant and hence also the 'cone' of nilpotent 
matrices A/2 = {X G M2(Z) | X 2 = 0} as described by Proposition 

The similarity class of the matrix X is [X] = {TXT~ l \ T G GL2(Z)}. 

Let E = ^ q ^ ; Notice that E l £ [E] using conjugation by ^ ^\ ; also 

—E £ [E] using conjugation by the matrix D = f ^ ^ 

Let N denote the set of non-negative integers. Let E\ = [XE], A G N; it 
is clear that E\ = X[E]. Also, it is an easy check to see that £\ and for 
A //x, are disjoint. 

Proposition 2.2. A matrix X £ A/2 is similar to XE, for a unique X G N. 
Thus A/2 = Uagn is the disjoint union of similarity classes. 

Proof. If any entry of the matrix X is zero, then it follows immediately 



from Proposition 2.1 that X is either the zero matrix, a multiple of E, or a 



multiple of E l . Thus, in this case X is similar to either the zero matrix or 
E. 

Suppose now that none of the entries of the matrix X is zero; say X = 
x y 



. Since the determinant is zero, the non-zero rows are dependent 

z —x J 

over the rationals; furthermore, after factoring out the greatest common 
divisor, say A, of the entries of X, we find that there are relatively prime 
integers m and n so that mx = nz and my = — nx. 

Now, using the relatively prime condition as above, we obtain the divis- 
ibility: m\z, m\x, n\x, n\y. Rewriting then x = mnxi, y = ra/i, z = mz\\ 
canceling where possible then from the original equations, it follows that 
mxi = z\, —nx\ = y±; thus m\zi and n\y\. Hence, x = mn\,y = —n 2 X,z = 
m 2 X, AGZ. 

Thus the original matrix is a multiple of a matrix, which has a factoriza- 
tion, 



(ran -n 2 \ = ( n\ , _ 
ym 2 —mn J Im J ^ ' ' 



Consequently, the action of conjugation is just a mixture of the usual action 
on Z 2 as row and as column vectors: 



TXT' 1 =T[ n J (m -n)T-\ 



-l _^ I n 

The conjugation action of T = ( a ^ I gives 



TXT' 1 = ± 



' (an + bm)(cn + dm) —(an + bm) 2 

(cn + dm) 2 —(an + bm)(cn + dm) / 

Now it is easy to see, using the Euclidean algorithm, that for any relatively 
prime integers m,n there is an T E SL^Z), the elements of GL^Z) of 

determinant one, so that T ( ] = (J J. Hence using the conjugation of 



K m J \0y 

such a T on X as described above gives the matrix E. Thus any non-zero 
integer matrix X with relatively prime entries is similar to E. ■ 



3. Enumerating Pythagorean Triples 

Consider the matrices in M2(Z) which are similar to E; there are relatively 
prime integers m, n so that the matrix is 

(mn -n 2 \ _ 1 / C S - N\ 
\m 2 -mn)~2\S + N -C )' 

The integers S, C, N give a Pythagorean triple satisfying S 2 + C 2 = N 2 . 

We now have a method of generating Pythagorean triples. Consider the 
conjugation action of GL2(Z) on E and determine the orbit. Of course 
this is just the same as the set of cosets GL 2 (Z) by the stabilizer of E. 
The stabilizer in GL^Z) of the matrix E is easily seen to be the subgroup 
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generated by U and the scalar matrix —I. However, the situation is a bit 
subtle. These cosets parameterize the different matrices on the left side of 
the equation above, but not the different Pythagorean triples, which are 
obtained from the right hand side. For primitive (m, n relatively prime) 
Pythagorean triples S, C, N we only want to count the really different 
solutions; fix m to be even and n to be odd and assume both m and n 
are positive; this specifies C as even and positive, N is positive, and S is 
odd. To enumerate the distinct Pythagorean triples as described above, we 
shall use matrices of determinant 1 with even 2 — 1 entry, m; it follows 
that the 1 — 1 entry n is odd since the matrix T is invertible; hence also 

9 ) . Since the 

m —mn J 

effect of the matrix —I, the negative identity, (on triples S, C, N) is trivial, 
we mod out that action. So we have either both m and n are positive or 
opposite signs; since both situations occur the counting the cosets will double 
count the Pythagorean triples. We want to use the cosets to enumerate the 
Pythagorean triples in a 1-1 fashion; and so we look for a way avoid this 
duplication. We shall develop the group theory a bit more to aid in this. 

The homomorphism r : PSL^Z) — > PSL2(Z2), induced by reducing mod 
2, has kernel denoted, T(2). Since the image of r is the group PSL2(Z2) (it 
is a non-abelian group of order 6) the group T(2) is of index 6 in PSL2(Z) 
and normal. It is easy to see that T(2) is generated by (the images of) U 2 = 

J ^\ and L 2 = ^2 l) ' Moreover, these elements generate the group 

as a free product. This will be useful for the enumeration of Pythagorean 
triples. For any relatively prime integers n, m, one odd and one even, 
there is a matrix in PSL2(Z) with those as the entries of the first column; 
but multiplying by an appropriate coset representative, we have that these 
elements n, m (up to sign) are the first column of a matrix in T(2). 

Proposition 3.1. T(2) is the free product of the infinite cyclic subgroups 
generated by U 2 and V 2 . 

Proof. This result is easily proven using the structure of PSL2CZ) as 
the free product [jl], of the subgroups of orders 2 and 3 generated by A = 

° x ^ and B = Since L = AB, L 2 = ABAB, U = AB' 1 , 

U 2 = AB~ 1 AB~ 1 , then any alternating word in C/ ±2 or V ±2 is also alter- 
nating in A, B^ 1 and hence non-trivial. ■ 



The enumeration of triples is based on the enumeration of cosets of T(2) 
modulo the stabilizer of E, which is the subgroup generated by U 2 . Because 
of the free product structure the coset representatives are just the words in 
L 2 , U 2 ; since we are considering cosets of the subgroup generated by U 2 , 
the rightmost letter of a coset representative is L ±2 . We do a breadth first 
enumeration, adding {7 ±2 or L ±2 on the left of a string (alternating); for 



example, if the last group element on the left is L 2 we can add L 2 , U 2 or 

u- 2 . 

For any matrix M 6 M2(Z), let 5(M) be the matrix with the same off 
diagonal entries as M and with signs changed on the diagonal. We want 
to avoid counting both M and 5{M) in the enumeration of cosets of T(2) 
since they produce the same Pythagorean triple from the action on E. The 
conjugation action of the matrix D and multiplication by —I has the same 
effect as S, 8{M) = -(DMD~ V ) = DMD~ X (projectively). Thus 5 gives an 
outer automorphism of order 2 of T(2); conjugation by D. On the generators 
of T(2) the automorphism is: 

5(U 2 ) = U~ 2 ,5(L 2 ) = L~ 2 . 

We then can effectively describe the distinct coset representatives which 
will enumerate the distinct Pythagorean triples using the alternating prop- 
erty for the words of the free group T(2) in the generators L 2 and U 2 . 
Also, using 5 we see that to obtain distinct triples we can initialize with 
Cq = {L 2 }, Cq = {}, Uq = {} and define inductively, i > 0, 

£t +1 = {L 2 X\X€£+UUt}, 
Cr +1 = {L- 2 X\XeCT^jut}, 

u+ +1 = {u 2 x\xeutucf}, 

K w = {U- 2 X\X€UruCt}. 

Thus we have an effective enumeration of non-trivial primitive Pythagorean 
triples using the coset representatives in the disjoint union 

i>0 

The recurrence relation which counts the number of elements pk in V of 
level k is pk = 3pk-i, since we can add any one of 3 elements to the left of 
a given alternating string to keep it alternating. We use II to denote the 
triples obtained after conjugation of E by the elements of V . 

For each T 6 II, the entries determine a unique primitive Pythagorean 
triple of positive integers as follows: 

\S\ = |T 2) i + Ti )2 |,|C| = |2T 1;1 |,|iV| = |T 2;1 -T h2 \. 

Theorem 3.2. The primitive Pythagorean triples are in 1-1 correspondence 
with the set II of certain of the the coset representatives of the subgroup 
generated by U 2 in T(2). This set II can be enumerated as a union of triples 
of level k indexed by subsets Hk of size 3 k . 

To make the tree we connect the j th element of level i with the 3 elements 
of level 2+1 numbered 3(j — 1) + 1, 3(j — 1)+2, 3(j — 1)+3 corresponding to the 
3 possible elements used to make the next level by the conjugation action of 
three of L~ 2 ,L 2 , U~ 2 , U 2 ; this yields triples, up to sign, respectively, which 
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can further be exprssed in terms of the coordinates S = m 2 — n 2 ,C = 
2mn, N = m 2 + n 2 , 

r 0~\ 

L = [m — Amn + 3n , 2mn — An ,m — Amn + 5n J 
L + = [m 2 + 4mn + 3n 2 , 2mn + An 2 ,m 2 + 4mn + 5n 2 ] 

r O O O ~\ 

U = [— n + 4mn — 3m , 2mn — 4m , n — 4mn + 5m J 

r O O O ~\ 

U+ = [— n — Amn — 3m , 2mn + 4m , n + 4mn + 5m J. 

Amazingly, we only have to compute the three transformations L+, U_ , U+ 
when the coordinates S, C, N are all positive to get the next level in the tree. 
When S, C, N are all positive then the previous triangle is L_ . This is easy to 
see, since in this case m > n > 0, and therefore 2mn > 2mn — An 2 > —2mn 
so that the second coordinate of L_, 2mn — An 2 , is smaller in absolute value 
than C = 2mn. Similarly, the second coordinate of each of L + ,U_,U+ 
is greater than C. Furthermore, with variations on this kind of argument, 
we can modify these three transformations to obtain three new transforma- 
tions, so that all the coordinates are positive, and so preserve the positive 
quadrant of the cone x 2 + y 2 = z 2 . 



U- 



m 2 + 4mn + 3n 2 , 2mn + An 2 ,m 2 - 
S,-C,N]+2(N-S + C)[1,1,1] 



Amn + 5n 



n 



Amn + 3m 2 , Am 2 



2mn, n 2 — Amn + 5m 2 ] 
■S, C,N]+2(N + S-C)[l,l,l] 

0~\ 

n + Amn + 3m , 2mn + 4m , n + 4mn + 5m J 
S,-C, N] + 2(N + S + C) [1,1,1]. 



One can easily check that each of the differences between coordinates in 
[m 2 — n 2 , 2mn, m 2 + n ], namely, (m — n) 2 , 2n 2 , 2mn — m 2 + n 2 , persists to 
the next level (except for a sign change on the last one) for exactly one of 
U- , C+ , U+ . These three new points are so represented in the tree from left 
to right (see Figure |l|) but connected to the level i point. For example one 
can follow the sequence with difference 1 between the z and y coordinates, 
in the tree as the path 

[3, 4, 5], [5, 12, 13], [7,24, 25], [9, 40, 41],- •• . 

The difference 1 between the y and x coordinates however will alternate, 

[3, 4, 5], [21, 20, 29], [119, 120, 169], [697, 696, 985], 

In fact any difference at level % will continue up the levels of the tree in a 
similar fashion. The difference —7 occurs (alternating sign) for 

[15, 8, 17], [65, 72, 97], [403, 396, 565], 

One can ask for a characterization of all those integers which can be 
differences. This is just a question about the representations, for m even and 



n odd, of integers by the forms P = (m—n) 2 , Q = 2n 2 , R = {n-\-m) 2 —2m 2 or 
-R= [jn-nf-ln 2 . The first can take on any odd square value, the second, 
any twice an odd square, while the third case is more restrictive. The norm 
form X 2 — 2Y 2 = D can be solved for Y and X for exactly those integers 
D which are norms in the ring of integers for the field Q(\/2). However we 
require that Y is even, so this is more restrictive; since Y is even, consider 
x 2 — 8y 2 = D; then it also follows that any solution when D is odd has x 
even. Consider this equation modulo p for any prime divisor of D then it 
follows that 2 is a quadratic residue modulo p. Using quadratic reciprocity, 
it follows that every prime divisor of D must be ±1 mod 8. However, for 
a prime divisor of D, p = —1 mod 8, it follows that x 2 = p = —1 mod 8; 
which however is impossible. It follows that any odd prime factor of D 
having odd exponent must be congruent to 1 modulo 8. Since the ring of 
integers in Q(\/2) is a unique factorization domain it follows that for every 
p dividing D there is an element of norm p when p = 1 mod 8 (not just a 
multiple). Multiplying these together gives a solution for D. Furthermore, 
for the minimal positive representation of D by this form R we obtain an 
infinite path in the tree, since we can solve for m and n from the equations 
X = n + m,Y = m. 

Given the absolute values of x, y, solutions to x 2 — 2y 2 = R > 0, then we 
can solve for values of m, n yielding a Pythagorean triple. Choose n + m = 
x > and m = y < then n = x — m > 0. Hence the Pythagorean 
triple of positive integers [n 2 — m 2 , -2mn, m 2 + n 2 ] has first difference R = 
n 2 —m 2 —(—2mn). For example to find the first difference 17, one has the root 
Pythagorean triple coming from the solution to 5 2 — 2-2 2 = 17; giving rise to 
the solution, 5 + 2y/2, so m = —2, n + m = 5, yields the triple S = 45, C = 
28, N = 53; to solve for D = —17, we can multiply by the smallest unit of 
negative norm, 

-1+V% obtain the solution (-l + v / 2)(5 + 2 v / 2) = -l+3>/2 
so solving —R, we obtain n = 3, m — n = 1, which yields [7, 24, 25]. One can 
also use a method of Lagrange as described in Q, to solve x 2 — 8y 2 = p. For 
example when p = 89, one finds the smallest zq < 89/2, so that (z 2 — 8)/89, 
is an integer, say q\\ here zq = 39, q\ = 17; then solve for the smallest 
z\ < 17/2, with (zf — 8)/17 = c/2- In this case Z\ = 5, t/2 = 1 < and 
the algorithm terminates. It now follows that (39 - 2y/2)/(5 + 2^/2) has 
norm 89. This gives rise to the solution 11 + 4\/2, an element of norm 
89; the number 3 — 7\/2 has norm -89. The solution 11 + 4^2, yields 
m = -4, n + m = 11, so n = 15 yields the triple S = 209, C = 120, N = 241; 
from the solution 3 — 7\/2 to —R, we obtain n = 7, m — n = 3, which yields 
the triple [51, 140, 149]. The last difference is 9 which occurs in Figure |l| at 
[33,56,65]. As a check one may compute the three neighbors of this point 
yielding [209,120,241], [51,140,149], [275,252,373]. Thus [51,140,149] is 
the root of those points with first difference -89. 
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